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Abstract 

The external properties of the sun represent strong boundary conditions for 
possible models of the solar interior. Even the variation of the particle diffusion 
is restricted by the observed intensity of high energy solar neutrinos 


1. Introduction 

The possibilities for alternative solar models are limited by the external 
properties of the sun: The radial distribution of the mass density and of 
the total pressure are determined by the assumption of static equilibrium, 
by the multi-pole oscillation frequencies, the radius and the mass of the 
sun. The chemical composition of the sun is characterized on the one 
hand by the initial mass abundances of hydrogen, helium and heavier 
elements. On the other hand the central region of the sun has a deficiency 
of 4.6 x 10 55 protons and contains 1.16 x 10 55 additional ct-particles, 
which have been produced since the formation of the solar system. For 
a known distribution of hydrogen and helium the temperature and the 
radiation power L H h of the various proton-proton fusion processes can 
be evaluated. The contribution of the CNO cycles is proportional to the 
mass abundance 


X (> = X(C) + X(N) 
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of carbon plus nitrogen. This number is determined by the total radiation 
power L of the sun 

L hc = L - L hh <x X 6 , L = 3.846(1) x 10 26 W, 
where, here and later, the notation 

3.846(1) = 3.846 =t 0.001 

is used. 

In this way many parameters of a solar model are fixed already. Only the 
radial distribution for a given number of a-particles can be varied. Tem¬ 
perature and energy spectrum of the solar neutrinos depend sensitively 
on a change of the ct-particle distribution by diffusion. In this connection 
the neutrinos are important because from the solar radiation power one 
can deduce a mean neutrino flux of 6.54 X 10 14 m~ 2 s _1 Experimentally 
roughly one half of this amount can be detected. In the case of the chlor¬ 
ine and gallium experiments the reduction can be interpreted as a neu¬ 
trino oscillation. This means that a part of the electron neutrinos 
produced in the sun are transformed to muon and tauon neutrinos, which 
cannot induce an inverse beta decay in Cl- and Ga-nuclei on the earth. 
But in the Kamiokande experiment based on the electron-neutrino scat¬ 
tering there is no strong dependence on the special type of neutrino; in 
this experiment there is a preference to the high energy neutrinos from 
the B decay. The boron production is sensitive to the central temperature 
and thus sensitive to the helium distribution in the central region of the 
sun. In this case a reduced neutrino flux indicates a central temperature 
smaller than the corresponding value in the standard solar model without 
diffusion. In a previous contribution [1] temperature and diffusion have 
been estimated. In this paper the helium distribution function shall be 
discussed in a more systematic way. All contributions to the total pressure 
and all HH fussion processes shall be taken into account. Finally the pos¬ 
sible variations of a solar model receive a further restriction. Section 3 
shows the consequences for the diffusion of helium in the radial region, 
where the thermonuclear reactions are concentrated. In Section 4 the par¬ 
ticle distributions are varied and the results are summarized. A compar¬ 
ison of the expected neutrino fluxes with the experimental results favours 
a special distribution of the a-particles in the sun. 

2. The Radial Distribution of Helium 

The mass abundances X , Y and 1 — X— Y of hydrogen ^H), helium 
( He) and the other nuclides ( 2 H, 3 He, atomic number Z > 2) respec¬ 
tively, start with the initial values 

X 0 = 0.727, Y 0 = 0.256, 1 - X 0 - Y 0 = 0.017 (1) 
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The thermonuclear reactions in the sun led to the present mass abun¬ 
dances 

X = X 0 -Y Y=Y 0 + Y 1-X-Y = 0.017. (2) 

They are characterized by the mass abundance Y of 4 He, which has been 
produced during the history of the sun. The distribution of the heavier 
elements does not change by the HH and HC processes. The present par¬ 
ticle number densities n p of protons and n a of ct-particles can be expressed 
by the mass density p 

n p = 5.97538 x 10 26 kg ~ x pX, n a = 1.50456 x 10 26 kg“V^ (3) 

Since 4.64 Ga tr the sun is a star on the main sequence of the Hertzsprung- 
Russel diagram. As a time unit the tropical year 2000 is used 

la tr = 3.15569254 X 10 7 s, lGa tr = 10 9 a tr . 


It shall be assumed that the electromagnetic radiation power L,, increased 
from 72 to 100 per cent of its present value 


L, = 3.846 x 10 26 WI 0.72 + 0.28 


4.64Ga„ 


Then one gets a total number 

N a = J n a dV 


= 1.50456 x 10 26 kg“ 1 / pYdV — 1.156 x 10 ; 


"I- 


( 4 ) 


( 5 ) 


of additional a- particles, if the early production of thermonuclear energy 
is dominated by HH processes and HC processes can be neglected. It can 
be shown numerically that the specific radiation power is of the type 


dL r 

~dM 


= 1.82(7) X 10~ 3 ^ exp ( - 40 x 5/3 ), 


x = r/Rj, Rj- = 6.9626 x 10 8 m, 


( 6 ) 


where r is the distance from the solar center. jL,. means the radiation power 
at the distance r from the center. R s is the solar radius. The thermonuclear 
reactions are restricted to the central region of the sun (x < 0.36). 

It can be assumed that the mass abundance of the additional a-particles 
can be expressed by the function (6) 

Y = Yj (t3 ) exp ( - /?x 5 / 3 ) for x < 0.36 

Y = 0, for x > 0.36 


( 7 ) 
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where the constant Y t is fixed by the integral (5) 

/* 0.36 

Yi(/3) / p(x)[exp ( — /3 x 5 / 3 )]x^x = 18.108 kgm -3 (8) 

Jo 

In a solar model without diffusion (/? = 40), the central temperature has 
its maximum value. If diffusion is possible (/? < 40), then the central tem¬ 
perature becomes smaller. 

The choice of the mass density p is not arbitrary, but it is restricted by 
the fact that in an equilibrium state p and the gravitation pressure P are 
correlated with each other and with the velocity v s of sound 

r 2 dP = -G N M r pdr, M r = 4 tt f Q r p(s/R s )s 2 ds 

v s = (5P/3p) 1//2 , Gjv = 6.67259 x 10 -11 kg _1 m 3 s -2 , (9) 

where G N is Newton’s constant. Eq. (9) for v s holds in a totally ionized 
material. This velocity can be evaluated from the frequencies of solar mul¬ 
tipole vibrations. A mass density distribution, for which these values of v s 
can be reproduced, has been proposed by Stix [2]. For the central region 
of the sun (x < 0.36) the values of p , M r jM s , P and v s are given in Table 1, 
where M s = 1.9889 X 10 30 kg means the solar mass. The notation 

aHn = a X 10 n 


Table 1. 

Mass density p y mass M r , total pressure P and velocity v s 
region of the sun 

of sound in the central 

mi 

O/kgm 3 ) 

(M,./ Mj) 

(p/p») 

(Pj/kms -1 ) 

0.00 

1.530E5 

0 

2.372E16 

508.3 

0.02 

1.481E5 

8.531E-4 

2.310E16 

509.8 

0.04 

1.359E5 

6.472E-3 

2.142E16 

512.5 

0.06 

1.198E5 

2.019E-2 

1.905E16 

514.9 

0.08 

1.029E5 

4.344E-2 

1.642E16 

515.6 

0.10 

8.769E4 

7.626E-2 

1.380E16 

512.1 

0.12 

7.403E4 

0.1178 

1.136E16 

505.8 

0.14 

6.220E4 

0.1667 

9.197E15 

496.4 

0.16 

5.195E4 

0.2212 

7.334E15 

485.1 

0.18 

4.309E4 

0.2795 

5.751 El 5 

471.7 

0.20 

3.547E4 

0.3397 

4.455E15 

457.5 

0.22 

2.902E4 

0.4000 

3.415E15 

442.8 

0.24 

2.355E4 

0.4590 

2.602E15 

426.1 

0.26 

1.905E4 

0.5155 

1.970E15 

415.2 

0.28 

1.533E4 

0.5686 

1.484E15 

401.6 

0.30 

1.225E4 

0.6178 

1.U4E15 

389.3 

0.32 

9.750E3 

0.6626 

8.349E14 

377.8 

0.34 

7.756E3 

0.7030 

6.262E14 

366.8 

0.36 

6.176E3 

0.7391 

4.702E14 

356.2 
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Table 2. The mass abundances Xand Yof 'H and 4 He, the partial pressures P 7 , P/, and P e 
of radiation, ions and electrons, respectively, for (3 = 36 


mi 

X 

Y 

( Ty W 

(W 

(p,m 

0.00 

0.399 

0.584 

1.185E13 

1.016E16 

1.355E16 

0.02 

0.416 

0.567 

1.136E13 

9.957E15 

1.313E16 

0.04 

0.450 

0.533 

1.042E13 

9.347E15 

1.206E16 

0.06 

0.491 

0.492 

9.304E12 

8.436E15 

1.061E16 

0.08 

0.535 

0.448 

8.102E12 

7.370E15 

9.037El 5 

0.10 

0.576 

0.407 

6.683E12 

6.269El 5 

7.520E15 

0.12 

0.612 

0.371 

5.385E12 

5.215E15 

6.141E15 

0.14 

0.643 

0.340 

4.219E12 

4.255E15 

4.937E15 

0.16 

0.667 

0.316 

3.261E12 

3.415E15 

3.915E15 

0.18 

0.685 

0.298 

2.469E12 

2.691 El 5 

3.058E15 

0.20 

0.699 

0.284 

1.862E12 

2.092E15 

2.360E15 

0.22 

0.709 

0.274 

1.399E12 

1.609E15 

1.805E15 

0.24 

0.715 

0.268 

1.070E12 

1.229E15 

1.372E15 

0.26 

0.720 

0.263 

8.141 Ell 

9.320E14 

1.037El 5 

0.28 

0.723 

0.260 

6.362E11 

7.030E14 

7.799E14 

0.30 

0.724 

0.259 

4.835 Ell 

5.283El 4 

5.848E14 

0.32 

0.726 

0.257 

3.802E11 

3.965E14 

4.380E14 

0.34 

0.726 

0.257 

3.007 Ell 

2.977E14 

3.282E14 

0.36 

0.727 

0.256 

2.382E11 

2.237E14 

2.463E14 


is used. For a given value of (3 the mass abundances X and Y, the partial 
pressures P 7 , P h and P e of radiation, of baryons (nuclei, ions) and elec¬ 
trons, respectively, are determined for each distance r in the radial region 
under consideration. Table 2 shows the corresponding quantities for the 
parameter (3 = 36. In this region (x = 0. 0.36) by Pauli’s principle the 

pressure P e of the electrons increases by 3.9 . 0.9 per cent compared to 

a classical gas of electrons. 

Table 3 contains the temperatureT(in units of Megakelvin), the relative 
atomic mass and the ratio L r /L for (3 — 36. The temperature is fixed by 
the equilibrium condition that the sum 

P 7 Pb + P e = P (10) 

becomes equal to the gravitation pressure P. For the evaluation of the 
relative atomic mass [i it has been assumed that the gas is ionized totally 
For (3 = 36 a value X 6 = 0.0049 is needed to reproduce the (electromag¬ 
netic) radiation power L. If the neutrino emission is taken into account, 
then the total radiation power 

L tot = 3.940 x 10 26 W 

becomes by 2.45 per cent larger than the radiation power L. 


(ii) 
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Table 3. Temperature T, relative atomic mass (i, electromagnetic radiation power L r 
at the distance r from the center for /? = 36. Opacity k(P) for ft = 30, 36 and 40 in units 

of m 2 kg -1 


m.) 

(T/MK) 


(L,./L) 

k, (30) 

k (36) 

k, (40) 

0.00 

14.721 

0.8074 

0 

- 

- 

- 

0.02 

14.568 

0.7940 

0.0066 

- 

- 

- 

0.04 

14.258 

0.7684 

0.0531 

0.036 

0.078 

0.116 

0.06 

13.859 

0.7392 

0.1506 

0.045 

0.080 

0.108 

0.08 

13.388 

0.7110 

0.2889 

0.069 

0.098 

0.120 

0.10 

12.759 

0.6862 

0.4525 

0.095 

0.117 

0.132 

0.12 

12.088 

0.6657 

0.5961 

0.123 

0.137 

0.146 

0.14 

11.373 

0.6495 

0.7239 

0.149 

0.157 

0.160 

0.16 

10.664 

0.6371 

0.8234 

0.176 

0.178 

0.178 

0.18 

9.947 

0.6279 

0.8895 

0.202 

0.200 

0.197 

0.20 

9.270 

0.6214 

0.9378 

0.222 

0.217 

0.213 

0.22 

8.630 

0.6168 

0.9632 

0.236 

0.230 

0.226 

0.24 

8.070 

0.6137 

0.9811 

0.253 

0.246 

0.242 

0.26 

7.538 

0.6117 

0.9893 

0.266 

0.260 

0.256 

0.28 

7.087 

0.6103 

0.9950 

0.295 

0.290 

0.287 

0.30 

6.617 

0.6095 

0.9973 

0.319 

0.314 

0.312 

0.32 

6.231 

0.6090 

0.9989 

0.328 

0.324 

0.322 

0.34 

5.876 

0.6087 

0.9998 

0.361 

0.358 

0.356 

0.36 

5.544 

0.6085 

1 

- 

- 

- 


If the energy transport is due to radiation in thermal equilibrium, then 
the opacity k — k{(3 ) is given by the equation 


16tt 2 d 4 

U = --—ctsb r —T , 
3 k p or 


( 12 ) 


where ct S b means the Stefan-Boltzmann constant. Table 3 shows that the 
opacity does not depend on the parameter (3 for r > 0.16 R r For smaller 
distances r from the center 36) is not essentially smaller than «(40). 
Therefore, the radiation equilibrium does not appear as a special problem 
in the region /3 = 36 ... 40. 


3. Diffusion of a-Particles 

The number n[J of additional a-particles, which have been produced 
in a spherical volume of radius r, is given by the expression 

N^ } = 1.156 x 10 55 L r /L. 


(13) 
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The present number of additional a-particles in the same volume is equal 
to the integral 

Nor = 1.50456 X 10 26 kg _1 x 4t r [ Y (s/R s )ps 2 ds. (14) 

Jo 

Both numbers are shown in table 4 for (3 = 36. Here the considerations are 
restricted to the central region of the sun because 98 % of the radiation 
power and of the a-particle production arise from r < 0.24 R s . The differ¬ 
ence between the numbers (13) and (14) can be expressed by the particle 
flux :j a (;r/R. r , t ) of a-particles 

N^. - Nar = 47rr 2 / j a (r/R s ,t)dt, t\ = 4.64Ga tr 
Jo 

j a (r/R n t) = (0.72 + 0.28///, )j ai {r/R J ), (15) 


where(r/R s ) is the present flux of a-particles in radial direction 

N ( ar - Nar = 3.44*/,^,(r/R,). (16) 


Table 4 shows these values of the a-particle flux for (3 = 36. They may be 
compared with theoretical estimates of the diffusion current density j d 
evaluated from the change AA of the mean free path A in radial direction 


= — <t q . = 2.18 6: 


2a, 


Toe 


\ J M a vl 


A = 


1 


Up (7p + (y 


’ Va * M a ksT '■ 


1/2 


(17) 


Here a p and cr a , are the scattering cross sections of ap and aa scattering 
orthogonal to the radial direction, respectively, if the target nuclei are 
screened by a gas of free electrons, ayis the fine-structure constant. M a , 
means the mass and P a , the mean velocity of the a-particles. n e , n p and n a 


Table 4. Numbers and N ar of u-particles produced and now present, respectively, 
within a sphere of radius r. Present current density j a \ and estimated diffusion flux jj of 

a-particles for /3 = 36 


ms) 

N (0) 

a r 

Nor 

(y a1 /m- 2 s- 1 ) 

(+/m 2 s ') 

0.04 

6.139E53 

5.713E53 

3.475E19 

2.426E19 

0.08 

3.338E54 

3.100E54 

4.838E19 

2.608E19 

0.12 

6.888E54 

6.420E54 

4.237E19 

1.845E19 

0.16 

9.515E54 

9.034E54 

2.449E19 

1.228E19 

0.20 

1.084E55 

1.053E55 

9.902E18 

8.393E18 

0.24 

1.134E55 

1.121E55 

2.918E18 

5.965E18 
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are the particle number densities of electrons, protons and a-particles, 
respectively. Contributions of the order /M a v a ) have been neg¬ 

lected in the expression (17) for the cross sections. For the diffusion cur¬ 
rent density one finds 



For /? = 36 both current densities, j a i and1 jj, are in the same order of mag¬ 
nitude. Thus the numbers in Table 4 suggest a solar model with j3 — 36. 

4. Conclusion 

Possible variations of the solar model are restricted strongly by the exter¬ 
nal properties of the sun (mass, radiation power, multipole vibration fre¬ 
quencies, total number N a of ct-particles produced within the sun). 
Therefore, only the parameter (3 , which characterizes the present distri¬ 
bution of the a-particles, can be changed. The consequences of such a 
variation are shown in Table 5. For /5 = 30 . .36 the present current den¬ 
sity j a i is in the same order of magnitude as the theoretically estimated 
diffusion flux jj. The value /? = 40 however, can be excluded, because in 
this case the current densityis too small; it even becomes negative for 
r = 0.08 Rj.. The radiation power L H h arising from various pp cycles 
remains nearly unchanged as a function of the parameter / 3 . Thus the mass 
abundance X 6 in the central region of the sun is fixed by the experimental 
radiation power L. An effective vector coupling constant G' v for the pro- 

Table 5. Ratio of the a-particle current densities, the contribution Thu of the pp-pro- 
cesses to the radiation power, the mass abundance X 6 of carbon and nitrogen, the central 
temperature^., the total neutrino flux (p Ui flux (p (B) of neutrinos from the 8 B-decay, the 
neutrino reaction rates (pa for the target nuclei 37 C1 and 7l Ga. Dependence on the 

parameter f3 


0 

Unl/j**) 

(L,i,,/W) 

x t 

Ct/MI<) 

30 

0.5 ..2.3 

3.634E26 

82E-4 

13.848 

36 

0.5 ..2.3 

3.628E26 

49E-4 

14.721 

40 

-0.1... 0.6 

3.617E26 

34E-4 

15.411 

0 

OA/itTV 1 ) 

W>(B)/m' 2 s _1 ) 

(0cr(Cl)/snu) 

(0<r(Ga)/snu) 

30 

6.5426E14 

0.8E10 

3.6 

135 

36 

6.5431E14 

2.2E10 

5.1 

139 

40 

6.5439E14 

4.5E10 

7.7 

148 

exp. 

- 

2.7(5)E10 

2.6(4) 

70(8) 
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cess p(p> de + is) must be smaller than 1.03 Gy, otherwise one gets L HH > L 
even for X 6 = 0. For G' v = Gy and 0 = 36 ... 40 the corresponding 
values of X 6 cannot be excluded. For /3 = 30 the abundance X 6 becomes 
too high. Thus only /? = 36 can be accepted. 

This means that the central temperature becomes T c = 14.7 MK rather 
than 15.4 MK (0 = 40, standard solar model without diffusion). The total 
radiation power 

Ltot = 3.94 x 10 26 W 

includes the neutrino radiation too. L tot and the total flux 

4> v = 6.54 x 10 14 ^/m 2 s 

of solar neutrinos, both are nearly independent of the parameter /3. The 
special neutrino flux <fi(B ) arising from boron decay has been evaluated 
under the assumption 

m a t ; 6 

For (3 = 36 this flux is consistent with the experimental values (Kamio- 
kande collaboration [3]). For (3 = 30 the neutrino flux is too small; for 
0 = 40 it would be too large. As the Kamiokande experiment is based 
on the elastic ise scattering, which should be roughly independent of 
the special neutrino generation, a parameter 0 = 36 is favoured. 

The reaction rates 4>cr for the chlorine and the gallium experiment are 
given in solar nutrino units (snu = 10 - 36 s -1 per target nucleus). For 0 = 36 
both reaction rates (f)cr(C 1) and (f)cr(G a) reach one half of the experimental 
values (references [4] and [5]). The reduction of the reaction rates is an 
indication for finite rest masses of the neutrinos: Either the neutrinos 
have a finite magnetic moment, which reduces the energy and the reac¬ 
tion rates of them; or by neutrino oscillations the electron neutrinos are 
partially transformed to muon and tauon neutrinos, which cannot induce 
inverse beta decays of Cl and Ga isotopes. The Kamiokande experiment 
however, can be reproduced quantitatively by a solar model with a diffu¬ 
sion parameter 0 = 36. 
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